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ABSTRACT 


A calibration technique is proposed that will allow the calibration of certain 
angular measurement devices without requiring the use of an absolute standard. 
The technique assumes that the device to be calibrated has deterministic bias errors. 
A comparison device must be available that meets the same requirements. The two 
devices are compared; one device is then rotated with respect to the other and a 
second comparison is performed. If the data are reduced using the technique 
described below, the individual errors of the two devices can be determined. 


INTRODUCTION 

All normal calibration techniques, whether for length, voltage, pressure, resis- 
tance, etc., involve the comparison of the test device against a measurement system 
or standard of known accuracy. Standard practice requires that a measurement 
device must be compared to a calibrated standard having an accuracy at least four 
times better than the test device (ref. 1). This requirement means that ultimately the 
standard must be calibrated by another laboratory with higher accuracy devices until 
a national or international standards laboratory is reached. 

This is not always possible and it is often desirable to be able to determine the 
errors in a device when an adequate standard is not available. A technique has been 
developed that allows the calibration of angular measurement devices that have 
certain types of errors. If a second similar device is available for comparison, the 
errors in the two devices can be determined without previous knowledge of the 
errors in either device. If certain conditions are met, the technique described below 
shows how to extract individual errors in the two devices. 

A technique has been developed in which two angular measurement devices such 
as protractors, inclinometers, encoders, resolvers, and Hirth couplings (ref. 2, 3) can 
be compared and the individual errors extracted. The requirements are that both 
units be stable, the errors in each must be periodic, and a means must be available to 
compare the two devices. The last requirement may be simple such as the case when 
two encoders are coupled together, and the comparison involves simply setting the 
shafts at different angles and reading the outputs of the devices. In other cases, such 
as the comparison of two Hirth couplings, an external device such as an autocollima- 
tor or an accelerometer may be required to monitor the difference between the 
devices. 


ANALYSIS 

Assume two angular measurement devices are to be calibrated and that the errors 
in both are repeatable and periodic. This means the error in the two devices can be 
represented by: 


ej = Ajsin(0) + A2sin(20) + AjsinpQ) + . . . 
-I- B 1 cos(0) + B 2 cos(20) -l- B 3 cos(30) + . . . 
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( 2 ) 


and, 

e 2 = A 1 ’sin(9) + A 2 ’sin(20) + A 3 ’sin(30) + . . . 

+ B 1 ’cos(0) + B 2 ’cos(20) + B3’cos(30) + . . . 

Assume a method is available to measure the deviation between the two devices. 
The devices are oriented so that each starts at a known angle and a calibration is 
performed showing the differences between these devices. The set of deviation data 
obtained at various angles can be converted to the following approximation: 

* ' = (Al ' Al,) Si ” (8) + Si ” (2e) + (A3 ' A3,) Si " (39) + ' ' ' O) 

+ (Bj-Bj’) cos(0) + (B 2 -B 2 ) cos(20) + (B 3 -B 3 ’) cos(30) + . . . 

If one of the devices is rotated with respect to the other by a known amount $, 
and a second comparison calibration is performed, then the new deviation can be 
characterized by: 

S 2 = A 1 sin(0)-A 1 ’sin(e+$) + A 2 sin(20) - A 2 ’sin(2(0 +$)) 

(4) 

+ A 3 sin(30) - A 3 ’ sin(3(0 + $)) + . . . 

+ Bj cos(0) -B 1 ’cos(0 + $) + B 2 cos( 20) - B 2 ’cos(2(0 +$>)) 

+ B 3 cos(30)-B 3 ’cos(3(0 +$)) + . . . 

It can be shown that equation (4) can be reduced to: 

S 2 = sin(0) [Aj - cos($) + B x ’ sin($) ] 

+ sin(20) [A 2 - A 2 ’ cos(24>) + B 2 ’ sin(24>)] 

+ sin(30) [A 3 - A 3 ’ cos(3$) + B 3 ’sin(3$)] + . . . 

+ cos(0) [Bj - Bj’ cos(4>) - Aj’ sin($)] 

+ cos(20) [B 2 - B 2 ’ cos(2<t>) - A 2 ’ sin (2$)] 

+ cos(30) [B 3 - B 3 ’ cos(34>) - A 3 ’ sin(3$)] + . . . 


Depending on the amount of offset, 4, certain values in this series do not 
produce independent data. For example, if $ =90 degrees, then the 4th, 8th, 12th, 
etc., terms of the series produce the identical information in both sets of data so that 
the coefficients of these terms cannot be determined. The simplest method to 
estimate the number of harmonics, N, required to characterize the errors is to 
examine the initial set of deviations using a Fourier transform, determine the num- 
ber of harmonics needed to closely fit the data, and then chose the value of 4 by: 


$ = 360/(N + 1) degrees 


( 6 ) 
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The technique works best when there are large changes in the patterns from the 
two sets of measured data. The maximum change occurs when the value of $> is half 
of the spacing between the most pronounced cycle of error. Even though selecting a 
high vame for N allows the periodic data to be more precisely fit, care must be taken 
not to select N too large. If N is high then $ is small and the shifts in the comparison 
patterns are slight and the error caused by small random variations becomes large. 
In the calibration of angular devices N = 3 has been shown to work well. The 
approximation to the error sets is good, and setting $ = 90 degrees caused cons ider- 
able differences to appear between the two sets of data. 

If we choose $ = 90 and limit our approximation to the first three terms, the 
equation (5) can be simplified to: 


S2 ' (Al + Bl>) Sin(9) + <A2 + A2> Sin(29) + <AJ ' B3 ° Sin(39) (7) 

+ (B x - Aj’) cos(9) + (B 2 + B 2 ’)cos(20) + (B 3 + A 3 ’)cos(30) 

If the coefficients of equations (3) and (7) can be determined then the individual 
values of the original error terms can be calculated. The coefficients of the sets of 
measured deviation data can be determined by Fourier analysis, least squares regres- 
sion, or numerical approximation to continuous Fourier coefficients. Once these 
coefficients are known, there are 12 equations and 12 unknowns which will permit 
solving for the coefficients that characterize the errors of the two individual 
measurement devices. 


SAMPLE CALIBRATION 

A comparison was performed using two Hirth couplings, each having a manufac- 
turer’s accuracy specification of less than 1 arc second, as the test devices and a 
precision servo accelerometer as an indicator of deviation. The data are shown in 
table 1. The second set of data, taken with device 2 shifted 90 degrees, appears in 
column 3 of table 1. 

The initial set of data was approximated by a discrete Fourier transform algo- 
rithm. The first three terms of the series fit the data with a standard deviation of .06 
seconds (see figure 1). The second series was also fit and the standard deviation was 
also .06 seconds (see figure 2). 

The following equations represent the fitted curves. Note that the zero frequency 
component is included here although it has no meaning since the zero on the divid- 
ers is arbitrary. 


S 1 = .2889 + (-.1556) sin(e) + (1.3389) sin(2e) + (-.0111) sin(3e) 
+ (-.0889) cos(e) + (-.2833) cos(2e) + (.1500) cos(30) 

S 2 — -.1769 + (-.2895) sin(0) + (-.2929) sin(2o) + (-.1476) sin(30) 
+ (.1799) cos(0) + (.0019) cos(20) + (-.0091) cos(30) 
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We now have 12 equations and 12 unknowns that can be solved to find the equa- 
tions for the error in the two devices under test: 

A r A{ = -.1556 A 2 -A 2 ’ = 1.3389 A 3 -A 3 ’= -.0111 

Bj-B^ = -.0889 B 2 -B 2 = -.2883 B 3 -B 3 ’ = .1500 

Aj+ Bj’ = -.2895 A 2 + A 2 = -.2929 A 3 -B 3 = -.1476 

B r A{ = .1799 B 2 +B 2 ’= .0019 B 3 + A 3 ’ = -.0091 

This can be solved by matrix manipulation or by simple algebra to determine the 
following: 


Aj = -.357 

A{ = -.201 

A 2 = .523 

A 2 ’ = -.816 

A 3 - -.159 

A 3 ’ = -.148 

B x * -.021 

= .067 

B 2 = -.141 

B 2 ’ = .143 

B 3 = -139 

B 3 ’ = -.011 


Thus the error in device #1 is approximately: 

ei = -.357 sin(0) + .523 sin(2e) - .159 sin(3e) 

-.021 cos(8) - .141 cos(20) + .139 cos(3e) 

The approximate error in device 2 is: 

e 2 = -.201 sin(9) - .816 sin(2e) - .148 sin(3e) 

+ .067 cos(0) + .143 cos(2e) - .011 cos(3e) 

Figure 3 shows the calculated errors in the two test devices. 

PROOF OF CONCEPT 

A method was devised to demonstrate the ability of this technique to extract 
individual errors. The method involves the error extraction of two devices with 
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relatively large errors and the subsequent calibration of these devices against cali- 
brated standards. This demonstrates the ability of the technique to separate errors 
and provides calibration curves that can be compared directly with the calibration 
curves determined by normal means. The actual values of the errors are large but 
the technique is based on the shape of the curves so that magnitude is not critical. 
The ability to extract degree errors to a tenth of a degree would imply the ability to 
extract arc second errors to a tenth of an arc second. 

Two generalized "protractors" (radial spoke patterns with lines about every 10 
degrees) were drawn on paper (figures 4 and 5) with intentional errors built in. The 
errors in each were generated as 3-harmonic sinusoidal series with about 3 degrees 
jeak amplitude but were distinctly different in the two. The protractors were over- 
aid and aligned at zero degrees. The differences were estimated by eye using a 
inear scale near the circumference, and are shown in column 2 of table 2. One 
protractor was then rotated 90 degrees (to the ninth line) with respect to the other 
and a second set of deviations was recorded (table 2, column 3). The estimated 
error in making these comparisons was about .25 degrees. The deviation plots are 
shown in figures 6 and 7. 

The data were reduced by the technique described in this paper and the individ- 
ual errors extracted. The two protractors were each calibrated against a normal 
protractor to determine their true errors. The accuracy of this comparison was 
about .25 degrees. The derived curves are compared with the true values in figures 8 
and 9. 

As can be clearly seen the correction curves determined in the two methods 
correlate very well. The standard deviation of the differences between the errors 
using the direct techniques and the new technique was .31 degrees. This proves the 
concept described in this paper. 


CONCLUSIONS 

A technique has been developed that will allow certain devices to be calibrated 
against each other without the need for a standard. If the errors in each are repeat- 
able, periodic, and able to be characterized by a reasonable number of sinusoidal 
harmonics, the individual errors in the two devices can be determined. 

The exact effects of scatter and drift in the test devices and the effect of errors in 
the comparison measurement have not been thoroughly investigated. This remains 
as possible future work in this area. Other possible future work involves using this 
method to intercompare three high-accuracy devices. This would allow the extrac- 
tion of the data on a common device by comparison to two other devices. This 
would further verify the method and show how much error should be expected when 
using this technique when comparing very high accuracy devices. 
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Reference angle 




■Em SBH 


0 

0.00 

0.00 

10 

0.44 

- 0.26 

20 

0.95 

- 0.47 

30 

1.29 

- 0.62 

40 

1.28 

- 0.55 

50 

1.27 

- 0.55 

60 

1.24 

- 0.60 

70 

1.08 

- 0.51 

80 

0.82 

- 0.47 

90 

0.47 

- 0.32 

100 

0.05 

- 0.20 ' 

110 

- 0.40 

- 0.25 

120 

- 0.72 


130 

- 0.91 

■■■■S3 

140 

- 0.96 

- 0.32 

150 

- 0.96 

- 0.32 

160 

- 0.84 

- 0.37 

170 

- 0.50 

- 0.32 

180 

- 0.09 

- 0.35 

190 

0.39 

- 0.34 

200 

0.85 

- 0.35 

210 

1.45 

- 0.32 

220 

1.89 

- 0.34 

230 

1.96 

- 0.30 

240 

2.00 

- 0.19 

250 

1.66 

- 0.15 

260 

1.29 

- 0.14 

270 

0.73 

- 0.07 

280 

0.10 

0.07 

290 

- 0.44 

0.18 

300 

- 0.72 

0.34 

310 

- 0.92 

0.63 

320 

- 1.10 

0.64 

330 

- 1.02 

0.59 

340 

- 0.82 

0.42 

350 

- 0.43 

0.14 


Table 1. Deviation between Hirth couplers. 
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Figure 1 . Comparison of unshifted indexers. 
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Figure 2. Comparison of shifted indexers. 
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Figure 3. Extracted errors in individual indexers. 
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Figure 6. Deviation between protractors (unshifted). 
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Figure 7. Deviation between protractors (#1 shifted 90 degrees). 
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Figure 8. Results of calibrating protractor #1 . 
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